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Abstract Polymer rings in solution are either permanently entangled or not. 
Permanent topological restrictions give rise to additional entropic interactions apart 
from the ones arising due to mere chain flexibility or excluded volume. Conversely, 
entangled polymer rings systems may be formed by closing randomly entangled flex- 
ible linear chains. The dependance of linking numbers between randomly entangled 
rings on the chain length, more specifically the second topological moment {n^)-, 
i.e. the average squared linking number, may be determined. In this paper, an 
approach recently discussed in mathematical physics and called abelian BP theory, 
is presented which allows to express the linking constraint in its simplest form, the 
Gauss integral, in terms of two gauge fields. The model of Brereton and Shah for 
a single ring entangled with many other surrounding rings is rederived. The latter 
model is finally used to calculate the second topological moment, in agreement with 
a recent result by Ferrari, Kleinert, and Lazzizzera obtained by using n-component 
0^ theory in the limit n — 0. 
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1 Introduction 



The topology of knots and links has attracted the interest of mathematicians and 
physicists alike (see for a general introduction e.g. |I|, 0, ^, ^, ^, ||). In polymer 
physics, the role of topological constraints has been discussed since its beginnings 
0, §. On the level of a self-avoiding walk description for the polymer conformation, 
the Gauss invariant (defined below) was discussed by Edwards who also considered 
the role of higher order link invariants (HOLC). In the abelian Chern-Simons 
theory (however not under the present term which is taken from the authors of 
||1CI|| ) was used to express the self-linking number of a single ring. In his discussion 
of HOLC, Edwards pointed to the role of a 3-vertex connecting 3 gauge fields, a 
remark anticipating the non-abelian Chern-Simons theory which was used by Wit- 



ten [11] some 20 years later to obtain polynomial invariants for knots and links. 



These algebraic invariants which are usually defined by mathematical knot theo- 
rists in terms of recursion relations (skein relations) with respect to the crossings 
appearing in the 2-dimensional projections of knots and links have also been applied 
to polymers (for this important approach which will not explored here, the reader 
is referred to 0). The conjecture by Edwards has been validated by the work of 
Cotta-Ramusino, Guadagnini, et al. (see 0, |^ and references therein) and by Ax- 



elrod et al., and by Alvarez and Labastida (see ||T3[ and references therein), who 
show that the Gauss invariant and all HOLC appear in a perturbation series of the 
non-abelian Chern-Simons theory when the vacuum expectation value of products 
of Wilson loops is evaluated. This result suggests to treat the Gauss invariant as an 
approximate measure of the true topology. 

The concept of the abelian variant of Chern-Simons field theory was reconsidered 
by many authors, among them |]12|, |I4|, ^ |1^, [1^, also for many rings systems, as 
it arises naturally whenever the Gaussian linking number is enforced as a constraint 
in the partition function. Partial reviews of the subject are H, IS, H]. The non-abelian 



variant has been discussed only recently by some authors [0, ^, 15], however its 
implementation in terms of a well motivated physical model for entangled polymers 
is not so clear. 

Whenever the abelian Chern-Simons theory is used in the context of the Gauss 
integral, an important problem arises in terms of self-linking integrals, where the 
double line integration is not taken with respect to two different ring contours but 
with respect to a single ring contour. These integrals attain a meaning only when 
using the so-called vertical framing, i.e. by replacing the single ring by a band. For 
details and the relevant literature we refer the reader to ]0]. 

Following a suggestion of Thompson (see e.g. the review ]|2^), an alternative 
topological quantum field theory, termed BF theory ]^, is used in the present 
paper. In contrast to Chern-Simons theory, it avoids the appearance of self-linking 
integrals, and was introduced, at least to the knowledge of the author, to the subject 
of polymer topology in ]jT5|. Recently, BF theory was rediscovered in |1^, |19| which 
is also concerned about avoiding self-linking integrals. In agreement with existing 
literature on topological quantum field theory as done in j^, the term BF theory 
will be used to underline the difference to Chern-Simons theory explained above. 

In the present paper, the direct link model for a system of Up rings of is 
used. It singles out a single ring, say Ci and enforces a linking number n of this 
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test ring with all other rings Cp G {C2, ■ ■ ■ , C*np} in the system. The Gauss invariant 
expressing the linking number as a double line integral is then represented in terms 
of an abelian BF theory. Next, the conformational degrees of freedom of all rings 
but Ci are summed over, giving rise to a "mass term" to be added to the BF 
theory. Then, one gauge field is integrated over, and the field theory of Brereton 
and Shah is obtained, expressed in terms of one remaining gauge field and the 
conformational field r^(s), where s parametrizes the polymer contour. This theory is 
then used to evaluate what has recently been called the second topological moment 
p!8[] , (n^). Essentially, the purpose of this paper is two-fold: one, to inform polymer 
physicists of a new approach (BF theory) to tackle entangled polymers, and two, to 
present a very simple derivation of the second topological moment in dense systems. 

The outline of the paper is as follows. In section 2, the direct link model is 
presented and expressed as an abelian BF theory. The model of Brereton and Shah 
is derived. In section 3, the second topological moment is evaluated. Finally, in 
the concluding section the relation to similar work is discussed and a simple scaling 
argument is evoked to interpret the result. 



2 The direct link model as a BF theory 

Let the conformation of the single test ring and all other rings be given by their 
coordinates {r^(s)} and {r^(s)} where P = 2, . . . ,np respectively. Now let n be the 
linking number of the ring Ci with all other rings Cf3 

Up 

n = Y,HCi,Cp). (1) 

where 

HC^,c^) = ^{ i dv-Adv^. (2) 



JcJc, |r--r/3|3 

is the Gaussian linking number. As discussed in the literature the double line 

integral on the r.h.s. is a topological invariant, i.e. is invariant with respect to all 
deformations of the rings given by their coordinates {r"} and {r^} which do not 
alter the topological state. However, this invariant is not unique. A counterexample 
is the Whitehead link |]T| which represents a linked state between two rings and has 
linking number (^{Ca-,Cp) = 0, the same value which results for a pair of unlinked 
rings. However, the Whitehead is maintained by self-interactions of one ring with 
itself. For random walk rings, these self-interactions effectively disappear, and thus 
the Whitehead link is unknotted to a pair of unlinked rings. 

The partition function of the system with a given linking number n reads as: 

Z{n) = ZoUn,J2HCi,Cp))j (3) 

The Kronecker delta appearing in Eq.(^ may be expressed in terms of the integral 

Sin, f) = r ^e^^-^^f (4) 



3 



giving rise to a topological "charge" g conjugate to n. One obtains the corresponding 
conjugate partition function Z{g): 



Z{g) = (expizgJ2HC,,C^)\j (5) 

The parameter g may therefore be considered as chemical potential for linking num- 
bers. In what follows, the conformational model for the rings is restricted to closed 
random walk chains. The excluded volume effect is neglected. As discussed in [|15|, 
this simplification restores the uniqueness of the Gauss invariant for links. 

The argument inside the average in Eq.(^) may now be expressed (in the spirit of 
a Hubbard- Stratonovich transformation) in terms of an abelian BF theory as follows 

w 



Z{g) = Af jvA j VBexp (^i j BAdA^ 



expitg f dr^ ■ A)) (exp WS-j ) (6) 

V JCi J I {ri(s)} \ V ^ // {r'3(s)} 

The constant is a normalization. The wedge product A is an exterior product 
which generates the algebra of differential forms over R^, f2*(M^) [^. In this lan- 
guage A and B are 1-forms, and the invariant expression B A dA is a 3-form. In 
terms of local coordinates it reads as e\^vB\d^Aydx^ A dx^ A dx^ . Therefore, no 
volume element appears in the integration. Line integrations are specified explicitly. 
The notation given above is widely used in the field theory literature , and is used 
here to make manifest the coordinate invariance of the field theory. Later on, when 
evaluating 2-point functions, we return to explicit coordinate representations of the 
fields. The variable j(x) is the tangent vector density, simply called the "current", 
of the chains /? = L.n^ 

\{^) = Y,f dsv\s)b{^-v\s)). (7) 

The conformational averages with respect to the test ring C\ and with respect 
to the other rings factorize for random walk chains. The conformation of the sur- 
rounding rings are represented in terms of the current j. In fact, they have fused 
to a single effective ring. As a reminder, the appearance of self-linking numbers has 
been avoided. 

Next, the conformational coordinates for the surrounding rings C\ for /3 7^ 1 are 
summed over. Then, the functional integral with respect to the gauge field B is 
performed. First, the average with respect to the ring chains Cp is carried out. For 
random walk rings one obtains: 

<^exp [i^B-^^ ^ = exp (^-^ j d^ j cix'B^(x)B,(x')(j^(x)j,(x'))) . (8) 

On the r.h.s. of this equation the chain index has been omitted. The average of 
the correlation function for the currents at points x and x' is carried out in the 
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limit of very long rings. Let us denote the number of segments per ring by and 
the average segment length by /. Neglecting terms arising due to the closure 



constraint |^ for the surrounding rings, the current- current correlation function is 



given by the following expression (the same approximation is made in 

(j,(x)j.(x')) = V^^(x-x) (9) 

where p = {npN)/V is the average segment density. Eq.(P) is an approximation 
valid for sufficiently concentrated systems. The question may be raised for these 
systems whether the Gauss invariant is still valid. Certainly the probability for 
more complicated links, such as the Borromean rings (a 3 ring link which falls apart 
when cutting one ring), which are only detected by HOLC, is higher but remains low 
compared to the Gauss invariant measuring the pairwise entanglement of rings. In 
order to determine the second topological moment given in the next section, these 
considerations are irrelevant as the Gauss invariant is simply an observable of the 
system (see below). 

Let G = pl'^/'i, then the partition function reads as 

Z{g) = Af jvA jvBexp(^i ^ B A c/A - ^ / ^ ' 

(expL (f dv'-A)) , (10) 

\ V JCi / I {ri(s)} 

where M is an adjusted normalization factor. Now the integration with respect to 
the gauge fields B may be carried out, which yields the effective model of Brereton 
und Shah p2i: 



Z{g) =^ J VAexp(^-^ J{V AA)'^ (exp(^tg £ dr' ■ A^"^ (11) 

Leaving aside for the moment the average with respect to the ring Ci and using the 
language of quantum field theory, one is left with a so-called Wilson loop that is 
averaged with respect to a euclidean Yang- Mills theory in 2 + 1 dimensions. Wilson 
loops are considered in quantum field theory to study the confinement problem. An 
analogy to this problem has been used in to study the collapse transition of 
randomly entangled polymer rings in 2D (for earlier studies see |2^, |2 



3 The second topological moment 

The second topological moment (n^) for a single test ring entangled with many other 
rings given the restriction that all conformational averages are taken with respect 
to random walks may now be determined rather simply. Following [|18], |19[ (n^) may 
be determined from: 

{n') = -^,Z{g)\,=o (12) 
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More generally Z{g) is the generating function for all topological moments. It is 
formulated for the present purpose as follows: 



zig) 



exp ( tg (p dr"'A 

Ca / / {A},{r"(s)} 



(13) 



The normalization factor has been absorbed into the average with respect to A. To 
separate the gauge fields and the conformational coordinates, the line integral in 
Eg .(1131) is expressed as follows: 



I dv-K = [ A^{k) (f ds r;.(s)e'''"(^) 

Jc Jk Jc 



(14) 



The summation convention is understood. The abbreviation /^represents J d^k/{2 
In order to perform the functional integral with respect to A, a gauge has to be used, 
which in the present case is V ■ A = 0. In this case the 2-point correlation function 
for A in Fourier space is given by: 



(A^(k)A^(q))|A} = 5(k + q)— ( 5^^ - k^K 



G 



(15) 



where = A;^/|k| are the components of a unit vector. In order to perform the 
conformational average, the only expression needed is 

$M.(si, S2; k) = (r^(si)r^(s2)e*'^-W^^)-'-(^^))) (16) 

This average has been evaluated for random walks rings ETf and reads as follows: 



<l>^^(si,S2;k) 



-V [Sis, - .2) - - 1 + ^Wsi -S2\{1- 



X exp 



Sl - -52 1 
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(17) 



The /c-dependent terms in the first factor on the RHS of this equation do not give 
any contributions due to gauge invariance. This is most obvious in the Landau 
gauge, k ■ A(k) = 0. 

Now, the second topological moment is given by (using a summation convention 
for repeated indices): 



in') = - 




A,{k)AM) d, dsi 6 ds2r^{si)r,{s2)e'''''^''^e'^<''^ ) (18) 



'k^q J J / {A},{r(.)} 

Using the correlators (0) and (|T^, one obtains the expression: 

1 



X exp 



G I — j) dsi J) ds2 

/2/.2 , |Si - S2 

Si - S2 1 



/2 / 1 

-^{S^,. - 1) U(si - S2) - — 
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(19) 



The last equation gives rise to 2 terms. The first stems from the delta function 
6{si — S2) and is readily evaluated. The second needs further analysis: if Si = S2, 
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the exponential equals 1, and one needs to evaluate an integral in k space which 
depends on a cutoff A ~ if si ^ S2, the integration with respect to Si, S2 
gives a non-trivial contribution to the integration in k-space, the details of which 
are discussed in the appendix. The result is when inserting for G 

, , NIA VN lA , 

In^) = pl^ — / 20 

^ ^ ^ ' 97r2 ^/2(3^ 97r2 ' ^ ' 

where A = C,^^ is an inverse cutoff length. As the result has been obtained in the 
limit ^ 1, the last term which is independant of A^ may be dropped, so the final 
result in the large A^ limit reads as: 

Apart from the numerical prefactors, the result is identical to the one derived from 



an n-component field theory in the limit n — > in [18|. The leading scaling term 



in Eq. (PT|) was also found by Brereton and Shah pH] who calculated the generating 



function Z{g) directly in the conformational space without the use of gauge fields, 
using a pre-averaging procedure as a further approximation. Their result contains 
however the square root of the ratio //^ contrast to the present result which is a 
minor difference and might be due to the approximation used in their work. 



4 Discussion 

The topological constraint of fixed linking number of a test ring entangled with Up 
surrounding rings has been implemented within the framework of a random walk 
model for the polymer conformation, by using the simplest link invariant, the Gauss 
integral. The latter one has been reformulated in terms of a topological quantum 
field theory, the so-called abelian BF theory. Two results have been obtained: first, 
the model of Brereton and Shah has been derived, second, the second topological 



moment in the limit of large segment numbers A^ ^ 1 has been calculated, in 
agreement with the previous result by Ferrari, Kleinert, and Lazzizzera [|l^, [T9 |. 



The same result may also be obtained from directly averaging the squared Gauss 
integral with respect to a random walk conformation, a work which will be presented 
elsewhere. The gauge field approach, however, is conceptually very appeahng as it 
separates topological interactions from conformational entropy. 

The present method avoids the complicated n-component field theory used in 



18]. The latter approach is, however, better suited to treat the full problem in- 



cluding excluded volume interactions. On the other hand, the Gauss invariant is 



ambiguous for self-avoiding walk rings. Considering work by Moroz and Kamien [30 
it remains an open problem whether a topological field coupled to the n-component 
field for the polymer conformation changes the random walk result significantly. The 
authors consider self- avoiding walks with writhe. In fact, a chemical potential for 
writhe is introduced which gives the coupling constant for the interaction between a 
topologocial gauge field and conformational n-component field, in a very similar way 
as in [|TB|. The scaling behavior of the radius of gyration and the first two moments 
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of the writhe are calculated. As to the first question, the effect of writhe is found 
to be irrelevant (to one-loop) , i.e. the self-avoiding walk fix point corresponding 
to the exponent v ~ 0.588 remains unchanged. Essentially this result is due to the 
fact that a topological field theory lacks a scale. Concerning moments of the writhe 
(which is similar to calculating moments of n as done above), no scaling dependence 
on the number of polymer segments is found which is a surprising result. A cleari- 
fication of this issue is certainly necessary. Let us note that Kholodenko and Vilgis 
3T|| calculated the writhe of semiflexible polymers which is found to scale as A^^/^. 

Finally let us interpret the result Eq.(^) in terms of a simple picture. Obvious 
the second topological moment basically scales as (n^) ~ pl^N as A^ becomes large. 
A simple argument, also discussed in [^], gives a similar result. Let the average 
segment density p inside the test ring be given by the number of rings crossing the 
interior of the test ring, which is equal to the average linking number per ring n, 
times the number of segments A^, divided by the volume of the test ring B? . Then 
the density reads as 

Assuming that the segment density in the interior of the test ring is equal to the 
segment density elsewhere in the system, and employing i?^ ~ /^A^^/^, i.e. using the 
Gaussian result for the ring volume, one immediately obtains 



n = pl^VN (23) 

Both n and \/Jn^ scale as a/ZV, so the theoretical calculation presented above repro- 
duces the characteristic scaling with chain length derived from the simple argument. 
The difference between density-dependent prefactors appears to be superficial. In 
fact in the case of dense melts where p ~ and which is the regime where the 
above argument applies, the prefactors coincide, and n = \/Jn?)- 
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Appendix 



In the following, the calculation of the r.h.s. of Eq.(|T9D 



(n^) = -G [ (f dsi (f ds2 



X exp 



\Sl - S2 1 



(24) 



6 ' 'V N 

is given in more detail. The delta function in the integrand gives rise to the integral 



3 Ak' 3 (27r)3,o 
GAT/^A (25) 
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where the k integration has been carried out up to the cutoff parameter A ~ ^ ^. 
The term in the integrand on the r.h.s. of Eq.(E4|) is given by 



21,2 



I'k 



6 



\si - S2I 1 



(26) 



Concerning the integration with respect to Si and S2 in the last equation, two cases 
need to be considered: Si = S2 and si 7^ S2. In the first case, the exponential 
function gives 1 and the integration gives a factor of A^, the result being 



'2a 



A 

3^ 



(27) 



The case Si 7^ S2 is slightly more involved. First, let us integrate with respect to k, 
which gives: 



'2b 



-G 



iV(2^ 



dSi (p dS2V2TT 



\si - S2I 1 



\Sl - S2\ 

N 



-1/2 



f28) 



Now the domain of integration with respect to Si and S2 may be restricted to two 
times one half the original one giving: 



l2b 



N 











^(^1-^2)11 



-1/2 



(29) 



A shift of variables u = {si — S2)/N, v = {si + S2)/N is easily performed and leads 
to 

2-a rp 



1 2b 



du 



dv 



3 A^(27r)3 

From now on the integration is elementary, and one obtains 

1 



-Nuil-u) 
3 ^ ' 



-1/2 



'26 



-GIN^/^ 



(30) 



(31) 



Adding Ji, /2a, and J26 and inserting for G = pl^/d for d = 3 gives the r.h.s. of 
Eq.(|19D in the main text. 
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